Localization Length Exponent, Critical Conductance Distribution and Multifractality 
in Hierarchical Network Models for the Quantum Hall Effect 
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We study hierarchical network models which have recently been introduced to approximate the 
Chalker-Coddington model for the integer quantum Hall effect (A.G. Galstyan and M.E. Raikh, 
PRB 56 1422 (1997); Arovas et al., PRB 56, 4751 (1997)). The hierarchical structure is due to a 
recursive method starting from a finite elementary cell. The localization-delocalization transition 
occurring in these models is displayed in the flow of the conductance distribution under increasing 
system size. We numerically determine this flow, calculate the critical conductance distribution, the 
critical exponent of the localization length, and the multifractal exponents of critical eigenstates. 

PACS numbers: 73.23.-b Mesoscopic systems; 73.40.Hm Quantum Hall effect (integer and frac- 
tional); 64.60.Ak Renormalization-group, fractal, and percolation studies of phase transitions 



I. INTRODUCTION 

Localization brought about by quantum interference in 
disordered media |u| is an ongoing research topic. One 
reason for this lies in the complex structure of the associ- 
ated localization-delocalization (LD) transition, a quan- 
tum phase transition which has no simple mean field de- 
scription H . Although many aspects of the LD transition 
are understood and quantitative results for critical expo- 
nents have been obtained by extensive numerical studies 
(see e.g. ||) there is no complete theory available yet. 
Best understood are quasi-one-dimensional (ID) systems, 
where the notion of a transfer-matrix led to strong an- 
alytic methods (see e.g. Hfl). Non-perturbative renor- 
malization techniques have also been developed on the 
basis of quasi- ID geometries. However, no such tech- 
nique proved to be successful directly for higher dimen- 
sional systems. There have been early attempts in that 
direction [^J?]] based on real-space renormalization within 
a Hamiltonian description of the systems dynamics. A 
major obstacle against their success was the increase 
in the number of relevant matrix elements of the effec- 
tive Hamiltonian after a few renormalization-group (RG) 
steps — unless the system was in a regime of strong local- 
ization. Only recently two groups have independently in- 
troduced a novel real-space renormalization technique for 
the LD problem |8|,[| . It was realized that a local scatter- 
ing theoretical formulation (see e.g. flf]|| ) of the systems 
dynamics gives rise to a straightforward real-space RG 
which is not restricted to quasi- ID. In particular, they 
considered the Chalker-Coddington (CC) model pd]|, a 
two-dimensional network of nodes and links where waves 
can propagate on the links and are scattered by chiral 
scatterers placed on the nodes. This model is very con- 
venient for testing the RG ideas since the scatterers are 
described by simple 2x2 random S-matrices. It has been 
established pT]-p~3[| that the CC model contains interest- 
ing LD physics. The chiral structure is meant to capture 



the essential features of the motion of 2D electrons in a 
stro ng p erpendicular magnetic field and a random poten- 
tial |11|] . Such electron systems exhibit the integer quan- 
tum Hall effect and are denoted as quantum Hall systems 
flilj . The quantum Hall effect is essentially caused by a 
LD transition Jl{| . In a localized state the electron is con- 
fined to some finite area within the bulk of the system 
while the electron's wave function is extended over the 
entire system in a delocalized state. However, states car- 
rying a unidirectional quantized edge-current exist also in 
the localized phase of a quantum Hall system. The link 
amplitudes of the CC model can be identified with local 
wave functions of a quantum Hall system. On increasing 
the average scattering strength (T) from to 1 one finds 
localization and chiral edge states (carrying a quantized 
Hall current) in almost any case. However, as soon as 
(T) is close to 0.5 extended bulk states appear, and e.g. 
the Hall conductance jumps to another quantized plateau 
value. It is expected that, in the thermodynamic limit, 
the transition occurs at precisely (T) — T* = 0.5. At 
the plateau to plateau transition the localization length 
£ of a critical state diverges, £ oc | (T) — T*\~ v , where the 
critical exponent v = 2.35 ± 0.03 [Q. Furthermore, it 
is known that the states are multifractal as soon as their 
localization length becomes larger than system size |l3| . 

It turned out, that the qualitative aspects of the 
LD transition could be described by a simple Migdal- 
Kadanoff type approach || leading to simple algebraic 
RG equations. Critical exponents, however, were found 
to be different from the commonly accepted values. 
Quantitatively more reliable results were obtained by l| 
when using a slightly more realistic approximation to the 
CC model. In that case, however, the treatment needed 
some numerical assistance. Despite the success in de- 
scribing the LD transition the RG approach was based 
on uncontrolled approximations. 

It was therefore tempting to elaborate further on the 
RG ideas of MM and develop a general scheme which en- 
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ables to improve quantitative results in a systematic man- 
ner. Advantageous for that purpose is the interpretation 
of the RG in terms of hierarchical networks embedded in 
the CC model. The hierarchical models are constructed 
from elementary cells cut out of the CC model and are 
characterized by the number of nodes, 6 2 , of the elemen- 
tary cell and a number of active scatterers, V. For each 
model the flow of the conductance distribution, the criti- 
cal conductance distribution and the critical exponent of 
the localization length can be determined by iterative RG 
steps. Furthermore, a calculation of multifractal proper- 
ties of critical states in the hierarchical models allows to 
complete the description of LD transition characteristics. 
It turns out that our RG approach leads to reasonable 
values already for moderate sizes of elementary cells. 

We have organized the paper as follows. In Sec. 2 we 
introduce the hierarchical models and describe their con- 
nection to those considered in |HJ. In Sec. 3 the flow 
of the conductance distribution is analyzed. The critical 
exponent of the localization length and the critical con- 
ductance distributions are determined. The multifractal 
analysis of critical states in hierarchical models forms the 
content of Sec. 4, followed by our conclusions in Sec. 5. 

II. HIERARCHICAL NETWORK MODELS 

We recall the construction of the Chalker-Coddington 
(CC) network model. A set of elementary scatterers is 
put on the nodes of a regular 2D lattice. Unidirectional 
channels link the scatterers to each other as shown in 
Fig. 1. The elementary scatterers are described by a 
unitary 2x2 scattering matrix S, 

T = \t\ 2 = \t'\ 2 is the transmission to the left and 
R = \r\ 2 = \r'\ 2 = 1 — T is the transmission to the right. 
The phases on the links are taken as random variables 
uniformly distributed between and 2tt. Also the trans- 
mission ('scattering strength') T E [0, 1] can be taken 
from a certain distribution. 

To study the localization-delocalization behavior of the 
CC network model one can impose boundary conditions 
for a square network of linear size L such that only total 
transmission to the left T(L) or right R{L) = 1 - T(L) 
is possible. In an appropriate two-probe conductance 
measurement (see Fig. 2) the transmission T(L) yields 
the conductance G via the Landauer-Biittiker theory of 
coherent transport, 

G=^T. (2) 

Thus, an interesting object to study the LD transition 
is the distribution function Pl{T) and its flow with in- 
creasing L. In the thermodynamic limit, L — » oo, bulk 



localized phases are characterized by T = or T = 1. 
In both cases current can only flow along edge states. 
Thus, two stable (attractive) fixed point distributions 
Pj(T) = S(T) and P„(T) = S(T - 1) exist. The LD 
transition will be characterized by a third fixed point 
distribution, denoted as P*(T), which is unstable (re- 
pulsive). The unstable character means that P* is the 
limiting distribution, lirar,^,^ Pl(T) = P*(T), only if 
the set of its relevant parameters Xi(L) (e.g. average 
value, median, variance) is initially taken to the fixed 
point values X*. In any other case the distribution ap- 
proaches one of the two stable fixed point distributions, 
Pi n. For large enough initial sizes Lq it may turn out 
that one parameter X(Lq) is enough to fix the flow of the 
whole distribution. This situation is referred to as one- 
parameter-scaling regime. The corresponding /3-function 
(3(X) — dX/d\nL describes the scaling flow of X(L) and 
allows for calculating the critical exponent v of the cor- 
relation length (localization length) jL7|, 

£oc \X-X*\- V , (3) 

from 



x* 

Therefore, a possible way of determining the critical 
exponent v would be to set up a renormalization-group 
(RG) process for the flow of Pl(T). An exact process 
could be achieved by recording the values T(L) for arbi- 
trary system sizes L and for an arbitrary large number of 
realizations of the random scatterers. Such process could 
be carried out with the help of numerical calculations and 
is only limited by the computational capacity. In a more 
efficient way, however, one replaces the network model 
by some hierarchically constructed network that allows 
to perform the renormalization recursively. The hierar- 
chical structure requires some finite elementary cell which 
allows for calculating the corresponding transmission T 
analytically or by simple numeric routines. However, the 
replacement of the original network by a hierarchical one 
introduces some error which is not under control. To 
reach definite answers for critical exponents one has to 
repeat the RG process for larger and larger elementary 
cells such that eventually the original model can be recov- 
ered. This procedure is advantageous if (i) one is aiming 
at qualitative results and (ii) the results for the critical 
exponents rapidly converge. To get qualitative results 
about the LD transition the elementary cell should be 
as simple as possible, still containing essentials like basic 
symmetries. For example, the chiral structure and paths 
with loops in the CC network model should be main- 
tained in order to allow for quantum interference along 
oriented paths. This idea was used in j|] where a simple 
elementary cell with 4 active scatterers (those for which 
neither T = nor T = 1) was considered (see Fig. 6 
in 0). In that approximation the scattering turned out 
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to obey one-dimensional composition laws which helped 
to proceed analytically. Unfortunately, in this approxi- 
mation a discrete mirror symmetry (with respect to the 
central vertical axis) is broken. As a result, the criti- 
cal distribution P*(T) is not symmetric with respect to 
T = 0.5, and the critical average transmission is different 
from 0.5. 

The obvious choice for an elementary cell respecting 
also the mirror symmetry is shown in Fig. 3a. It con- 
sists of a 3 x 3 lattice containing 5 active scatterers and 4 
scatterers (in the corners) that fully transmit to the left 
(T = 1) or to the right (T = 0). Building a hierarchi- 
cal network from such elementary cells leads to a fractal 
network with fractal dimension D^ 3 / 5 ! = In 5/ In 3 rj 1.46 
(see Fig. 3a). It will be denoted as [3/5]-model. The 
corresponding RG process has been already carried out 
in Ref. S. In that work a different interpretation was 
adopted. As in the [3/5]-model an elementary cell of 
5 scattering units is replaced by one in each renormal- 
ization step (see Fig. 2 in ||). As shown in Fig. 3b the 
scaling factor is 2 instead of 3. In other words, the RG 
process of proceeds by constructing a new network 
from renormalized scatterers while our process considers 
a hierarchical network as a fractal subset of the origi- 
nal CC network. Apart from the different interpretation 
both processes are, of course, equivalent if the difference 
in scaling factors is taken into account. The simplified 
elementary cell, considered in Fig. 6 of ||] corresponds to 
the [3/5]-network when replacing the central scatterer by 
a fully reflecting/transmitting one (see Fig. 3 c) and will 
henceforth be denoted as [3/4]-model. Its fractal dimen- 
sion is L>[ 3 / 4 ! = In 4/ In 3 w 1.26. 

In the present work we extend the RG process to larger 
elementary cells and study also the multifractal proper- 
ties of critical wave functions. 

A generalization to larger elementary cells is straight- 
forward: each Chalker-Coddington network with b x b 
sites (b = 5, 7, 9 . . .) and boundary conditions of only 4 
links to the outer world can serve as an elementary cell. 
By recursion one gets a hierarchical infinite network (see 
Fig. 3a, c). Its fractal dimension is 



where the elementary cell contains V = (b 2 + l)/2 scat- 
terers. Note, that for b — > oo the original CC model is 
recovered and e.g. the fractal dimension approaches 2 for 
6> 1. 

In general, hierarchical networks can be characterized 
by an odd number b describing the size of its elemen- 
tary cell (b 2 ) and the number V of active scatterers. Ac- 
cordingly we denote the hierarchical networks as [6/V]- 
network. 

In addition to the hierarchical networks with V(b) = 
(b 2 + l)/2 which, for large b, reduce to the original 
CC model one can use [6/U]-elementary cells to con- 
struct larger elementary cells of size b' x b' with V(b') < 



(b' 2 + l)/2 active scatterers, in order to simplify the RG 
steps. We elaborated on a [9/37]-model constructed from 
[3/5]-elementary cells (displayed in Fig. 8). Note that 
the number of active scatterers (V(9) = 37) is less then 
(9 2 + l)/2 = 41 which would follow from the straightfor- 
ward construction for 6 = 9. 



III. FLOW OF THE CONDUCTANCE 
DISTRIBUTION 

After having outlined which hierarchical models will be 
considered in this work we are now going to describe how 
the RG for the conductance distribution was obtained in 
our numerical calculations. 

1. In the first step we algebraically solve the scattering 
problem for a particular elementary cell. We obtain the 
transmission probability T(L = b) as a function of the 
scattering matrix elements t = \t\e la , t' — \t\e la , r = 
\r\e l13 and r' — \r\e l13 of all active scatterers. Here the 
unitarity constraints 

e »(a+o') = _ e *(/9+/9') 

and 

|t| 2 + |r| 2 = l 

have to be obeyed. That leaves for T(b) an algebraic 
expression of V scattering strengths, T\^y, and 3V 
phases. 

2. In the second step we draw the scattering strengths 
T from an initial distribution, Po(T) peaked around a 
certain mean value (T) and random uncorrelated phases 
for a, a' and /?. For each realization we calculated the 
scattering strength T(L = b) according to the algebraic 
expression. 

3. After collecting a large number (« 2000) of scatter- 
ing strengths we calculated the corresponding histogram 
and normalized it. That yields an approximation for the 
distribution function PL=b{T) depending on the initial 
distribution Pq(T). We also determined the average value 
(T) (L = b) and other quantities like median or geometric 
mean of the distribution. 

4. We repeat the steps No. 2 and 3, but now draw- 
ing T from the previously obtained distribution Pi,=b(T). 
We did this for each active scatterer independently. As 
a result we got a novel distribution which should coin- 
cide with PL=b 2 (T) of the hierarchical [6/U]-model if the 
phases at each RG step can be taken as uncorrelated 
random phases. Indeed, the phases of t(L = b),r(L = 
b),t'(L = b) turned out to be uniformly distributed in 
the interval [0, 2ir[ and correlations were not observed. 
Therefore, the 4th step yields a reasonable approxima- 
tion for P L=b i(T) and, in particular, for (T) (L = b 2 ). 

5. Next, we repeat step No. 4 N times (N 20) which 
yields, for each initial distribution Pq(T), the flow of the 
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distribution function Pl (T) and of the mean (T) (L) as 
functions of InL = (N + 2) In b. 

6. We repeated the RG process for a number of differ- 
ent initial distributions Pq(T). 

An example for the flow of the distribution is displayed 
in Fig. 4 for a [3/4]-system with initial value (T) < T c ps 
0.6. It demonstrates that the distribution is attracted by 
the localization fixed point at T = 0. 

It generally turned out that after a few iterations the 
mean (T) was already a good scaling variable, i.e. once 
a certain mean value was reached, the flow of the dis- 
tribution after further iterations did no longer depend 
significantly on the initial distribution. In particular we 
found for each model a certain distribution which was 
independent of L, characterized by (T) (L) = T* but un- 
stable against changes (T) - T* ^ 0. Thus, we took 
the mean (T) to investigate the instability of the critical 
distribution. The critical exponent v of the localization 
length is determined as 

= (JV-JVo)lnb 

V bx{{T{N)-T*)/(T(N )-T*)}' {) 

This expression follows from Eq. (||) with X = (T) . The 
numbers No and N were chosen such that (T) was sig- 
nificantly repelled from T* without leaving the regime of 
power law repulsion. 

In Figs. 5-8 the elementary cells together with the 
critical distribution are shown. For the [3/5], [9/37]- 
models also the flow of the average transmission is dis- 
played. The critical distribution turned out to be broadly 
distributed over the interval [0, 1]. Those models respect- 
ing the mirror symmetry of its elementary cell showed 
critical distributions approximately symmetric with re- 
spect to 0.5 ps T*. The [3/4]-model breaks the mirror 
symmetry and T* ps 0.6, as previously found in ||. This 
model and the [3/5]-model discussed in are taken as 
a reference to judge the improvement in quantitative re- 
sults when going to larger elementary cells. The symmet- 
ric critical distributions for different models are hardly 
distinguishable when taking numerical uncertainties into 
account. The critical distribution is not uniform over 
the interval [0, 1] (as conjectured in [p§), but develops a 
unique shallow minimum at 0.5 and has maximum val- 
ues at the extreme positions T = 0, 1. This confirms the 
observations of || for the [3/5]-model and of |l^] for the 
original CC model. 

The values for the critical exponent v are collected 
in Tab. 1. One observes that with increasing elemen- 
tary cells the critical exponent v approaches the CC 
network value in a monotonic way from above and the 
[9/37] value, v = 2.75 ± 0.3, is already close to the value 
v = 2.43 ± 0.18 obtained for the original CC model 
ffjj. Note that the value v = 3.5 ± 0.3 for the [3/51- 
model is compatible with the result [y ss 2.4) of ||] 
when taking the difference in scaling factors into account 
(In 3/ In 2 ps 1.58). 



IV. MULTIFRACTAL ANALYSIS OF CRITICAL 
STATES 

A rather complete characterization of a LD transition 
phenomenon in terms of critical exponents is provided by 
the knowledge of the critical exponent v of the localiza- 
tion length and the multifractal exponent ao of critical 
states |2j. The exponent ao describes how the typical 
value Pt of the probability P to find an electron (being 
in that state) in a small box of linear size I scales with 
respect to A = l/L, 

Pt(A)ocA a °. (7) 

The fact that critical states have ao > D in a D- 
dimensional system shows that critical states are ex- 
tended but not totally 'space filling' (multifractal). In- 
deed, the typical value scales faster to zero as the average 
value ((P) oc L~ D ) and thus their ratio 

p := P t / (P) oc L-^-°^ (8) 

is a convenient order-parameter for the LD transition. It 
vanishes algebraically on approaching the LD transition 
point from the extended phase and stays zero in the local- 
ized phase ||] . Also in the quantum Hall system, where 
no extended phase exists, ao essentially determines the 
scaling of critical states, including energetic and spatial 
correlations. 

The multifractal exponent ao has been numerically de- 
termined in [^3| for the original CC model and found to 
be 2.27 ± 0.01. The method used in that work can be 
described briefly as follows: The network model with its 
scattering matrices gives rise to a discrete unitary time 
evolution operator U operating on bond amplitudes and 
mapping a state ip(t) at some time t to a novel amplitude 
configuration ip(t + l) at one instant of time later by tak- 
ing all of the scattering conditions into account. Energy 
eigenstates ip are solutions of 

Utp = if) (9) 

and thus correspond to the eigenstates of U with eigen- 
value 1. The energy of the underlying electron system 
enters here only parametrically through the scattering 
matrices and U will have an eigenstate to eigenvalue 1 
only for a discrete set of the parameter T . Instead of 
looking for this discrete set one can fix T and ask for the 
more general eigenvalue problem 

U{T)^ n {T) = e^V^T) . (10) 

ipn{T) are, for fixed T, a set of eigenstates at scattering 
strength T of slightly modified disorder configurations, 
each of the modifications being an overall shift in random 
phases on the links. 

The typical probability -Pt(A) can be identified as the 
geometric mean of local squared amplitudes integrated 
over boxes of size I, 
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^t(A) :=cxp^ln(j]H 2 J^ . (11) 

Here the average (. . .} can be performed over different 
realizations and/or different local boxes within one real- 
ization. 

We adopted this method for computing the critical ex- 
ponent oq of the hierarchical network models. The scat- 
tering strengths were fixed to the critical value T*. We 
could not rely on a recursive calculation method, but 
needed some finite system of linear size L > b where b is 
the linear size of the elementary cell. This was necessary 
to detect the strong fluctuations of local amplitudes re- 
sponsible for the multifractal behavior. Instead of going 
to very large sizes L we took moderate sizes (L « 50) but 
used several realizations to obtain ao via Eqs. (|l^J^)- An 
example for a critical wave function of a [3/5]-system is 
shown in Fig. 9. Strong spatial correlations, previously 
observed in critical states ]2(i[ ] , are even more pronounced 
on the fractal support. 

Since the hierarchical networks are defined on geomet- 
ric fractals of dimension £)[ & / y ] it is inappropriate to com- 
pare different values of ao directly. Rather one should 
address the exponent of the order-parameter p which is 
ao — D^ b / V \ These values are listed in Tab. 1. It can be 
seen that ao — D^ h / y ^ monotonically approaches the CC 
model value from above as D^ b / y ^ approaches d = 2 from 
below. Since large values for ao — D correspond to strong 
multifractal fluctuations of critical amplitudes this be- 
havior tells that the hierarchical models with sparsely dis- 
tributed active scatterers show more pronounced fluctua- 
tions than those with a denser set of active scatterers. For 
the [9/37]-model, however, we obtain a^—D = 0.33±0.05 
which is already close to the CC value, 0.27 ± 0.01. 

V. CONCLUSIONS 

We extended the con- 

cept of a real-space renormalization-group approach to 
the Chalker-Coddington model m a systematic man- 
ner. We constructed hierarchical network models based 
on elementary cells which were cut out of a Chalker- 
Coddington model. Each elementary cell is characterized 
by the number b 2 of nodes and the number V of active 
scatterers. The hierarchical models are fractal objects of 
dimension Z)[ f, / 1/ ] — In Vj In b embedded in the original 
Chalker-Coddington model. In the V =oc b 2 — > oo limit 
the models coincide with the Chalker-Coddington model. 
We performed renormalization-group calculations for the 
following [6/^]-models: [3/4], [3/5], [5/13] and [9/37]. 
From these calculations we obtained the critical conduc- 
tance distribution and the critical exponent v of the local- 
ization length. The results for the [3/4] and [3/5]-models 
are consistent with those obtained in |^|5). Our results 
for larger elementary cells show how fast the critical 
exponent approaches the value of the original Chalker- 
Coddington model. 



Let us briefly comment on the possibility to use our 
method for a precise determination of the critical expo- 
nent of the original Chalker-Coddington model. For that 
purpose one not only has to use larger elementary cells, 
but also has to increase the number of realizations to 
get smoother histograms. Since the average value of T 
turned out to be a good scaling variable the precision of 
the critical parameter for each hierarchical model can be 
improved by increasing the number of realizations in each 
iteration step. The critical parameters arc then more 
precise for the particular hierarchical model, however 
have uncontrolled deviations from that of the Chalker- 
Coddington model. Better approximations can only be 
gained by increasing the size of the elementary cell and 
by studying the convergence of critical exponents. The 
computational effort then drastically increases; it is ex- 
pected to be comparable to that of the well known finite 
size scaling method based on the transfer matrix tech- 
nique [[|. 

To complete the characterization of the localization- 
delocalization transitions in the hierarchical [6/V^]- 
models we calculated the multifractal exponent ao — 
jj[ b / v \ f critical states and found that it approaches 
the Chalker-Coddington value from above when go- 
ing to larger elementary cells. In summary, our work 
shows that a systematic renormalization-group treat- 
ment of networks of random quantum scatterers can not 
only reveal qualitative aspects of a possible localization- 
delocalization transition, but also provides quantitative 
results for critical exponents the accuracy of which is con- 
trolled by convergence under stepwise increase of the size 
of elementary cells. 

This work was supported by the Sonderforschungsbere- 
ich 341 of the Deutsche Forschungsgemeinschaft. We 
thank Boris Shapiro for stimulating us to this work and 
Daniel P. Arovas for helpful discussions. MJ thanks 
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FIGURE CAPTIONS 



Figure 1: The graphical representation of the Chalker 
Coddington network model. Wave amplitudes propagat- 
ing on links can be scattered either to the left or to the 
right by unitary scattering matrices situated at the nodes 
of the network. 

Figure 2: Schematic representation of a Chalker- 
Coddington network of size L x L in a two-probe con- 
ductance measurement. The current flows from the 
source contact (down left) through the network to the 
drain contact (up right). The conductance is given as 
G(L) = (eyh)T(L). 

Figure 3: a) The construction of the [3/5] -hierarchical 
network. A 5 x 5 cell is cut out of the Chalker-Coddington 
model. Reducing the number of links to the outer world 
down to 4 links leads to an elementary cell with 5 active 
scatterers. The hierarchical construction yields a fractal 
network as a subset of the Chalker-Coddington model. 
At each level of the hierarchy the system size increases by 
a factor of 3. b) The renormalization scheme of Ref. [|| : 5 
scatterers are replaced by one. In each step of the renor- 
malization the rescaling factor is 2. c) The hierarchical 
model used in Ref. ||. It can be viewed as a subset of 
the [3/5]-model with fully transmitting/reflecting central 
scatterers. 

Figure 4: Flow of the distribution of transmission 
strengths under renormalization (in arbitrary units) . The 
initial distribution (top) is peaked around a mean-value 
below the critical value T*. After two renormalization 



steps (middle & bottom) the distribution clearly ap- 
proaches the stable distribution of full reflection (T = 0) . 

Figure 5: The elementary cell of the [3/4]-model (top) 
and the corresponding critical distribution of scattering 
strengths (bottom). 

Figure 6: The elementary cell of the [3/5]-model (top 
left), the corresponding critical distribution of scattering 
strengths (top right) and the corresponding flow of av- 
erage transmission strengths T under increasing number 
(N) of renormalization steps (bottom). 

Figure 7: The elementary cell of the [5/13]-model (top) 
and the corresponding critical distribution of scattering 
strengths (bottom). 

Figure 8: The elementary cell of the [9/37] -model (top 
left), the corresponding critical distribution of scattering 
strengths (top right) and the corresponding flow of aver- 
age transmission strengths under increasing number (JV) 
of renormalization steps (bottom). 

Figure 9: Squared amplitudes of a critical wave function 
for the [3/5]-model and a system size of L = 50. 



TABLE CAPTIONS 



Table 1: Table of models [6/V], fractal dimensions D, 
critical exponents of the localization length v and of the 
order-parameter exponents ctQ — D. 
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Model [b/V] 


D 




a - D 


[3/4] 


ft 1.26 


6.6 ± 0.6 


0.6 ± 0.05 


[3/5] 


ft 1.46 


3.5 ± 0.3 


0.4 ± 0.05 


[5/13] 


ft 1.59 


2.9 ± 0.3 


0.36 db 0.05 


[9/37] 


ft 1.64 


2.75 ± 0.3 


0.33 db 0.05 



Table 1: 



